Microscopic origin of exchange bias in core/shell nanoparticles 
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We report the results of Monte Carlo simulations with the aim to clarify the microscopic origin 
of exchange bias in the magnetization hysteresis loops of a model of individual core/shell nanopar- 
ticles. Increase of the exchange coupling across the core/shell interface leads to an enhancement 
of exchange bias and to an increasing asymmetry between the two branches of the loops which is 
due to different reversal mechanisms. A detailed study of the magnetic order of the interfacial spins 
shows compelling evidence that the existence of a net magnetization due to uncompensated spins at 
the shell interface is responsible for both phenomena and allows to quantify the loop shifts directly 
in terms of microscopic parameters with striking agreement with the macroscopic observed values. 
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I, INTRODUCTION 

Proximity between a ferromagnetic (FM) and an an- 
tiferromagnetic (AFM) material leads to interesting ef- 
fects that result from the structural modification and 
competition of different magnetic orderings at the in- 
terface between them. In particular, the exchange cou- 
pling at a FM/AFM interface may induce unidirectional 
anisotropy in the FM below the Neel temperature of 
the AFM, causing a shift in the hysteresis loop, a phe- 
nomenon known as exchange bias (EB). Although the 
first observations of this phenomenon, dating back five 
decades agoi, were reported on oxidized nanoparticles, 
most of the subsequent studies have focused on layered 
FM/AFM structures^ because of their application in 
advanced magnetic devices^. However, in recent years, 
the study of EB in nanoparticles and nanostructures^ 
has gained renewed interest since it has been shown that 
control of the core/shell interactions or of the exchange 
coupling between the particle surface and the embed- 
ding matrix can increase the superparamagnetic limit for 
their use as magnetic recording medial. Several exper- 
iments on different nanoparticle systems with oxidized 
shells have studied the sizel£, temperatur o^ft and cool- 
ing field dependence of the EB fields, as well as training 
effects^. However, the interpretation of the results may 
be hindered by collective effects and interactions with the 
embedding matrix since, up to date, no EB experiment 
has been conducted on a single particle, which would al- 
low to confront the results with the existing models. 

Both nanoparticles and layered systems display com- 
mon phenomenolgy, although in the later case a wider 
range of experimenal techniques have been used, which 
have provided deeper knowledge on the microscopic 
mechanisms that are at the basis of the EB effect. Re- 
cently, spectroscopic techniques have provided insight 
on the structure and magnetic behavior of the inter- 
face spins at a microscopic level, demonstrating the cru- 
cial role played by uncompensated interfacial spins on 
EBi^. Thus, knowledge of the magnetic structure at 
the interface has become a subject of primary interest in 
understanding EB. At difference from layered systems, 



the interface of core/shell nanoparticles naturally incor- 
porates roughness and non-compensation of the magne- 
tization, two of the main ingredients for which differ- 
ent assumptions are adopted by the existing models for 
EB in filmsiSiifi. Some microscopic models for bilay- 
ers have undertaken calculations of EB fields under cer- 
tain assumptionsiii^, numerical studies based on a mean 
field approach-^ or Monte Carlo (MC) simulation a 2 ft' 21 ' 22 
making different assumptions about the interface. How- 
ever, only very recently, some works partially address- 
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FIG. 1: (Color online) Hysteresis loops for a particle with ra- 
dius R = 12 a obtained from a ZFC state and after FC down to 
T = 0.1 in a field /ifc for Jsh = — 0.5Jc and Ji n t = — (+)0.5 Jc 
in the left (right) column. Panels (a) display the total nor- 
malized magnetization component along the field direction. 
Panels (b) show the average magnetization projection of the 
core spins along the field axis. Panels (c) show the normalized 
contributions of the shell spins at core/shell interface to the 
total magnetization of the loop. 
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ing the EB phenomenology in nanostructures have been 
publishe d 2 '^ 24 . In this article, we will show, through MC 
simulations based on a simple model of one core/shell 
nanoparticle, how some of the EB phenomenology to EB 
is related to exchange coupling at the core/shell inter- 
face. Moreover, the direct inspection of the magnetic 
configurations along the hysteresis loops will allow us to 
provide a quantitative understanding of the macroscopic 
loop shifts in terms of microscopic parameters. 

II. MODEL 

The nanoparticles considered have spherical shape 
with total radius R = 12a (a is the unit cell size) and 
are made of a FM core surrounded by an AF shell of 
constant thickness Rsh = 3a with magnetic properties 
different from the core as well as from the spins at the 
interface between core and shell spins. Taking a = 0.3 
nm, such a particle corresponds to typical real dimensions 
R ~ 4 nm and Rsh — 1 nm and contains 5575 spins, of 
which 45 % are on the surface. 

The simulations are based on the following Hamilto- 
nian: 

£ £ • / « X>< (s?) 2 -££-3 , a) 

(id) * i 

where Si are classical Heisenberg spins of unit magnitude 
placed at the nodes of a sc lattice. The first term is the 
nearest-neighbour exchange energy, where the value of 
the exchange constants Jy depends on the spins belong- 
ing to different particle regions. At the core, Jy is FM 
and will be fixed to Jc — + 10 K. Spins at the shell have 
AF coupling as corresponding to oxides; here a reduced 
value of the coupling at the shell with respect to the core, 
Jsh = — 0.5 Jc, has been set so that the Neel temperature 
of the AF TV is lower than the Curie temperature of the 
FM. Since, in real samples, it is difficult to access mi- 
croscopic information about the coupling at the interface 
spins [defined as those on the core (shell) with at least 
one neighbor on the shell(core)], we have considered dif- 
ferent values and signs for the exchange coupling at the 
interface, Ji n t. 

The second term accounts for the local uniaxial 
anisotropy along the z-axis. The anisotropy constant 
at the core is fixed to fee = 1 while the value at the 
shell, ksh = 10, is enhanced with respect to fee due to 
the reduced coordination of the shell spins. Finally, the 
last term describes the Zeeman coupling to an external 
field H applied along the easy-axis direction, which in 
reduced units reads h — fiH/ks- To simulate the hys- 
teresis loops, we use the MC method with a Metropolis 
algorithm. As for the spin updates, we use a combination 
of the trial steps which has proved useful for Heisenberg 
with finite anisotropies as described elsewhere 2 "^. Our 
protocol to simulate exchange bias mimics the experi- 
mental one: we first cool the system from a high temper- 
ature To > Tn disordered phase in constant steps down 



to the measurement temperature T — 0.1 in the pres- 
ence of a magnetic field h-pc — 4 K applied along the 
easy-axis direction. Then, the hysteresis loop is recorded 
using as starting configuration the one obtained after the 
field cooling (FC) process. The loops are obtained by cy- 
cling the magnetic field from h = 4 K to h = —4 K in 
steps Sh = —0.1 K and the different quantites averaged 
during 200 MC steps per spin at every field. 



III. RESULTS 

Typical hysteresis loops are shown in Fig. Q]for two 
values of the interface coupling Jj n t/Jc = —0.5, +0.5. 
Compared to the loops obtained from a zero field cooled 
(ZFC) state, the loops obtained after FC are shifted to- 
wards negative field values and have slightly increased 
coercivity (see Figs. QJi), independenly of the sign of the 
interfacial exchange coupling. In order to gain further 
insight on the differences between both cases, we have 
also computed the field dependence of the contribution 
of interface spins belonging the shell, MgjJ*, to the total 
magnetization as displayed in the lowest panels of Fig. 
For negative (positive) interface coupling, the interfacial 
spins at the shell acquire a negative (positive) net mag- 
netization after FC, in both cases higher than the one 
attained after ZFC, although more pronounced for the 
negative coupling case. These observations reflect that, 
after the FC process, a fraction of the interfacial spins 
(~15 % of the interface spins at the shell) have been 
pinned along a direction compatible with the core/shell 
exchange interaction, as corroborated also by the verti- 
cal shifts in the Mg^* loops (to be commented below). 
This is no longer true for the ZFC case, for which a high 
fraction of interfacial spins follows the reversal of the FM 
core as reflected by the change in sign of Mg^* along the 
hysteresis loop. 

Interesting enough, we have also noticed an increas- 
ing asymmetry of the FC loops with increasing values of 
the interface coupling, as it is apparent when compar- 
ing the descreasing and increasing branches of the loop 
in the top panels Fig. QJi- The origin of this asymme- 
try, is more clearly understood by looking at the average 
absolute value of the magnetization projection along the 
field axis through the reversal process, = J+ ' %\-> 
as depicted in Fig. for the core spins. This quantity 
presents peaks centered around the coercive fields that 
indicate deviations of the core magnetization from the 
applied field direction. In the ZFC case, the peaks are 
centered at similar field values and they are quite narrow 
and almost symmetric around the minimum. However, 
for the FC loops, apart for the obvious shift of the peak 
positions, the decreasing branch peak is symmetric and 
narrow while the increasing branch one is deeper and 
asymmetric, enclosing bigger area under the loop curve. 

These observations also indicate that the loop asym- 
metry reflects different reversal mechanisms in both 
branches of the hysteresis loops. This can be corrobo- 
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FIG. 2: (Color online) Snapshots of the spin configurations of a midplane cross section of the particle parallel to the z axis 
taken at fields of the descending and ascending branches close to the coercive fields {h^, Kq) for the case Ji n t = — 0.5Jc shown 
in Fig. la. 



rated by direct inspection of the spin configurations along 
the loops, as presented in the main panel of Fig. [21 for 
Jint = — 0.5 Jo As it is evidenced by the sequence of 
snapshots, the reversal proceeds by quasi uniform rota- 
tion along the descending branch, while nucleation of 
reversed domains at the interface and their subsequent 
propagation through the core center is basically the re- 
versal process along the ascending branch. Similar asym- 
metry between the loop branches has been also observed 
experimentally in bilayers2£i2L2&. A detailed inspection 
of the configurations, also reveal the presence of spins 
at core/shell interface aligned perpendicular to the field 
direction for intermediate field values (see for example 
the snapshots for h — —2.8, 1.0 in Fig. 0, This obser- 
vation corroborates the interpretation of recent results 
of small-angle neutron scattering experiments on Fe ox- 
idized nanoparticles in which the anisotropy of the ob- 
tained spectra was attributed to the existence of a net 
magnetic component aligned perpendicularly to the field 
direction22i2Si. 

The microscopic origin for the different reversal mech- 
anisms can be further claryfied by looking at the behav- 
ior of the interface shell spins along the hysteresis loop 
(see Fig. QJ:). While in the descending branch there is 
a considerable amount of unpinned spins that are able 
to reverse following the core reversal, in the ascending 
branch, remains constant (for Ji n t < 0), an indica- 

tion that spins at the shell interface remain pinned, hin- 
dering uniform rotation of the core but acting as a seed 
for the nucleation of reversed domains. The changes in 
the magnetic order at the core/shell interface and the 



presence of domain walls during reversal can be traced 
by monitoring the value of the average sum of the pro- 
jection of the spin direction into the direction of the total 
magnetization vector along the hysteresis loops as 

1 N 

m P (h) = -^2S i (h)-M i (h) . (2) 

i=l 

This quantity should be close to 1 if the magnetization 
reversal proceeds by uniform rotation of the spins, since 
in this case the spins remain parallel to the global magne- 
tization direction. Deviations of m p (h) from 1 indicate 
the formation of non-uniform structures during the re- 
versal process. An example of the field variation of m p 
computed for all the spins in the particle and for the in- 
terfacial spins is shown in Fig. where we have plotted 
separately the contribution of the core spins. 

During the decreasing field branch of the loop, m p re- 
mains quite close to 1 for the core spins, except for mod- 
erate decrease down to 0.7 for values of h close to the 
coercive field at this branch, . The sharpness and 
symmetry of the peak around confirms that the re- 
versal proceeds by uniform rotation. In constrast, during 
the increasing field branch, an increasing strong deviation 
of m p from 1 starting from negative field values can be 
clearly observed, reaching its maximum value also near 
the coercive field of the increasing field branch, h^, where 
m p ~ 0. In this case, the observed peak asymmetry is 
indicative of the nucleation of the non-uniform domains 
observed in the snapshots of Fig. El These domains are 
formed at those points of the core interface with weaker 
values of the local exchange fields, as indicated by the 
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FIG. 3: (Color online) (a) Shows the field dependence of the 
average spin projection into the total magnetization direction 
m p (squares) and the contribution of the core spins (circles). 
In (b), the contribution of all the interface spins (squares) has 
been taken into account, while the contributions of core and 
shell spins at the interface are shown in circles and triangles, 
respectively. 



more pronounced departure from 1 of mj, nt (/i) (see Fig. 
Et>), than those corresponding to the total magnetiza- 
tion (see Fig. Ek). The variation of m* 11 * for interface 
shell spins during the decreasing branch indicates the ex- 
istence of a fraction of shell spins that reverse dragged 
by the spins at the core, while constancy of m p in the 
ascending branch is indicative of spins pinned during the 
core reversal. 

Finally, we have also studied the variation of the coer- 
cive field and the EB field [defined as he — (h^ — h^)/2 
and /i e b = (h,Q + h c ~)/2 respectively] with the interface 
exchange coupling Ji n t, presented in Fig. 01 for positive 
and negative J\ n t values. For both J\ nt ^ 0, the values 
of he and /i e b are very similar, and a decrease in he and 
an increase in /i b is observed, with a nearly linear de- 
pendence for values of |Ji n t| smaller than the exchange 
coupling at the shell Jsh — — 0.5Jc- With the increase of 
| Jint |, core spins become more coupled to the unpinned 
shell spins, therefore facilitating the magnetazation re- 
versal with the subsequent decrease in the coercivity. At 
the same time, increasing | Ji n t| while keeping the values 
of Jc, Jsh and h-pc constant, results in higher local ex- 
change fields created by the uncompensated spins at the 
interface, causing an increase of the loop shift. Let us no- 
tice also that the values of the coercive and exchange bias 
fields obtained are of the correct order of magnitude when 
expressed in real units. For example, for Jint/ Jq = —0.5, 
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FIG. 4: (Color online) Variation of the coercive field he and 
the exchange bias field /i e (, with the exchange coupling con- 
stant at the core/shell interface for Jint < (open circles) 
and Jint > (open squares). The exchange bias fields com- 
puted from Eq. 0] as described in the text are shown with 
filled symbols for Jint < (down triangles) and Jint > (up 
triangles). 



we obtain He ~ 0.22 T and H c \, rs 0.11 T, which are in 
agreement with typical values found in studies of oxidized 
nanoparticleaViVi 1 "! 11 ! 12 ! 13 . 

The proportionality of h e h to J\ n t should be taken as 
a hint for the microscopic origin of the loop shifts. As 
we have mentioned in previous paragraphs, the observed 
vertical displacements of the loop corresponding to the 
interface shell spins point to the existence of a net magne- 
tization at the core/shell interface due to uncompensated 
pinned spins at the shell interface^. If this is the case, 
the coercive fields after FC can be thought as the sum 
of the ZFC coercive field h c and the local field acting on 
the core spins due to the net interface magnetization of 
the shell spins, so that they may be computed as 

h c = h° c + J Int M± t (3) 

and, therefore, the exchange bias field can be written as 

h eh = Jmt(M+ t + M£ t )/2, (4) 

where M Int = sh} z i^i ^ s the net magnetization 

of the interfacial shell spins at the positive (negative) co- 
ercive fields h,Q , and is the number of nearest neighbors 
of spin i. 

The last expression establishes a connection between 
the macroscopic loop shift and microscopic quantities 
that, although may not be directly measured in an ex- 
periment, can be computed independently from the sim- 
ulation results. The values of h e b obtained by inserting 
the Mj^ t values extracted from the Fig. QJ; in Eq. 0]are 
represented as filled symbols in Fig. 01 where we can see 
that the agreement with the h ch values taken from the 
hysteresis loop shift is excellent within error bars. 
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IV. CONCLUSIONS 

In summary, we have presented simulations of a model 
for onr core/shell nanoparticle that have revealed an 
asymmetry in the hysteresis loop due to the different 
magnetization reversal mechanisms in the two branches. 
This has been shown to be related to the exchange cou- 
pling at the interface J\ n t , independently of its sign and 
of any assumption about the nature of the interface re- 
gion. The detailed analysis of the changes in the mag- 
netic order of the interfacial spins have also allowed us to 
demonstrate that macroscopic EB fields can be computed 



microscopically from the knowledge of interface magne- 
tizations at the coercive fields. 
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